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Cover:  Helicopter  electromagnetic  sounding  over 
sea  Ice.  Hp  represents  the  primary  magnetic 
field  produced  by  the  transmit  coll  antenna  Tx. 
Hs  represents  the  secondary  magnetic  field  re¬ 
sulting  from  eddy  currents  caused  by  the  primary 
magnetic  field  flowing  In  a  conductor,  e.g.  sea¬ 
water.  Hs  Is  sensed  by  the  receiver  coll  Rx.  The 
In-phase  and  quadrature  response  of  Hs  Is  used 
to  estimate  the  distance  from  the  suspended  an¬ 
tenna  bird  assembly  to  the  seawater.  The  laser 
measures  the  distance  from  the  bird  to  the  snow 
or  Ice  surface.  From  these  determinations  the 
snow  and  Ice  thickness  Is  estimated. 
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Airborne  Electromagnetic  Sounding  of  Sea  Ice  Thickness 
and  Sub-ice  Bathymetry 


AUSTIN  KOVACS,  NICHOLAS  C.  VALLEAU  AND  J.  SCOTT  HOLLADAY 


INTRODUCTION 


In  May  1985,  a  multifrequency,  multicoil-geometry  airborne  electromagnetic  (AEM)  sys¬ 
tem  was  used  for  the  first  time  for  estimating  sea  ice  thickness,  water  conductivity,  and  water 
depths  from  about  1  to  20  m  under  an  ice  cover.  The  AEM  system  was  basically  a  standard 
geophysical  exploration  device  used  for  recognizing  and  analyzing  changes  in  the  earth’s 
conductivity.  The  technology  makes  use  of  the  fact  that  penetration  of  electromagnetic 
waves  into  conductive  materials  is  frequency-dependent,  as  are  the  return  amplitude  and 
phase  spectra  measured  at  a  receiving  antenna  (e.g.  Grant  and  West  1965.  Keller  and  Frisch- 
knecht  1966,  Kaufman  and  Keller  1983).  It  is  an  airborne  technology  that  has  been  used  by 
industry  for  several  decades  for  detecting  and  mapping  conductive  mineral  deposits  (e.g. 
Pemberton  1962,  Ward  1970,  Fraser  1981). 

The  concept  of  using  this  technology  for  measuring  coastal  bathymetry  was  reviewed  by 
Morrison  and  Becker  (1982)  and  for  detecting  sea  ice  thickness  by  Becker  et  al.  (1983).  The 
feasibility  of  using  an  AEM  system  for  measuring  sea  ice  thickness  perhaps  originated  in 
1968  (Anon.)  but  was  not  pursued  beyond  an  analytical  study.  However,  this  analysis  veri¬ 
fied  that  both  thin  and  thick  sea  ice  could  be  measured,  though  the  paper  cautioned  that 
pressure  ridges  would  pose  a  problem  because  of  their  three-dimensional  nature.  The  study 
also  mentioned  the  need  for  a  system  having  high  resolution  and  a  very  rapid  response  time, 
since  many  pack  ice  features  (for  example,  leads  and  ridges)  will  be  flown  over  in  seconds  or 
fractions  of  a  second. 

This  paper  presents  results  of  our  1985  study  which  determined  the  feasibility  of  using  an 
AEM  system  for  profiling  sea  ice  thickness,  seawater  conductivity  and  under-ice  water 
depth.  The  field  program  included  AEM  measurements  over  a  second-year  sea  ice  floe  of 
known  thickness  as  well  as  first-year  sea  ice  of  several  thicknesses.  Also  sounded  were  four 
very  long  and  thick  second-year  sea  ice  floebergs  grounded  on  a  shoal. 


Airborne  electromagnetic  system 

The  AEM  system  used  in  this  study  had  four 
pairs  of  coils  (transmitter  7*  and  receiver  Rx).  The 
coils  behave  as  magnetic  dipoles  and,  along  with 
system  electronics,  allow  simultaneous  operation 
at  four  distinct  frequencies  in  two  coil  configura¬ 
tions,  horizontal  coplanar  and  vertical  coaxial 
(Fig.  1).  The  nominal  operating  frequencies  were 
530  and  16,290  Hz  (horizontal  coplanar)  and  930 
and  4,160  Hz  (vertical  coplanar).  A  fifth  frequen- 
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Figure  1.  Three  typical  AEM  coil  ar¬ 
rangements.  Arrows  indicate  the  mag¬ 
netic  dipole  orientation. 


Figure  2.  Helicopter  electromagnetic  sounding  over  sea  ice. 


cy  of  32,020  Hz  was  also  evaluated  by  re¬ 
placing  the  16,290-Hz  coils  in  the  field. 

The  transmitter  and  receiver  coils  were 
about  6.5  m  apart  inside  a  Kevlar  tube 
(called  a  “bird”)  about  7.5  m  long  and 
0.5  m  in  diameter.  The  bird  weighed  about 
200  kg  and  was  typically  flown  about  35 
m  above  the  ice  surface  (Fig.  2). 

In  principle,  the  transmitter  coil  pro¬ 
duces  a  primary  magnetic  field  Hpf  which 
causes  eddy  currents  to  flow  in  nearby 
conductors  (e.g.  seawater).  The  resulting 
secondary  H,  and  primary  magnetic  Helds 
are  sensed  by  the  receiver  coils  (Fig.  3). 

The  distance  to  and  conductivity  of  a  con-  Figure  3.  Magnetic  fields  associated  with  AEM 
ductive  medium  affect  the  mutual  coup-  sensing  using  a  horizontal  coplanar  (whaletail) 
ling  ratio  HJHp.  Through  the  use  of  coil  arrangement. 
bucking  coils  in  the  bird  and  system  elec¬ 
tronics,  the  primary  Held  at  the  receiver  is 

cancelled  out,  and  highly  precise  measurements  of  the  in-phase  (IP)  and  quadrature  ( Q ) 
components  of  the  secondary  magnetic  Held  are  made  and  recorded.  The  measurements  are 
expressed  in  parts  per  million  (ppm)  of  the  primary  Held  that  would  have  been  sensed  by  the 


receiver  if  bucking  were  not  performed.  These  data,  along  with  frequency,  Tx  and  Rx  coil 
spacing  and  orientation,  and  bird  pitch  and  roll  sensor  data,  are  then  used  to  calculate  an  ap¬ 
parent  conductivity  a  and  bird  height  above  the  conductive  surface. 

2 


Since  sea  ice  is  relatively  resistive  and  therefore  transparent  at  the  AEM  system’s  lower  fre¬ 
quencies,  the  system  senses  the  conductive  seawater  and  thus  determines  the  distance  from 
the  bird  to  the  sea  surface.  A  laser  profilometer  system  was  built  into  the  bird  and  used  to 
measure  the  distance  from  the  bird  to  the  ice  surface.  Subtracting  this  distance  from  the  dis¬ 
tance  to  the  seawater  surface  determined  by  the  AEM  system  gives  the  apparent  ice  thick¬ 
ness,  or  the  snow  and  ice  thickness  where  a  snow  cover  exists. 

A  radar  altimeter  was  installed  on  the  helicopter,  and  a  radar  elevation  display  was  mount¬ 
ed  in  front  of  the  pilot  to  help  him  maintain  an  altitude  of  about  65  m  above  the  ice  surface. 
Since  the  bird  is  suspended  30  m  below  the  helicopter,  the  bird  was  flown  about  35  m  above 
the  ice  relief. 

A  flight-path  video  recorder  system  was  also  installed  on  the  helicopter.  The  camera  was 
mounted  to  view  the  bird  in  flight  and  the  terrain  being  overflown.  This  information  was 
available  for  post-flight  review  and  for  real-time  viewing  in  the  helicopter. 

General  theory 

There  are  a  number  of  ways  to  represent  the  responses  of  a  stratified  earth,  consisting  of  n 
layers,  to  an  oscillating  magnetic  dipole  elevated  above  the  earth’s  surface.  For  the  example 
in  Figure  4,  Tx  is  the  magnetic  source  dipole  at  height  /t, ,  Rx  is  the  magnetic  sensor  at  height 
h2,  x  is  the  horizontal  distance  between  the  electrical  center  of  the  Tx  and  Rx  coils,  and 
a,  .  .  .  oa  are  layer  conductivities  (S/m)  of  thickness  f,  .  .  .  fn.  Maxwell’s  equations  describe 
the  behavior  of  the  electromagnetic  field  as  follows  (for  a  general  case): 


V  •  (E  =  Qf 

(1) 

V  •  nH  =  0 

(2) 

VxE  +  jo>iiH  =  0 

(3) 

VxH  —  jo)(E  =  J 

(4) 

where  V«  =  vector  divergence 

Vx  =  curl  of  the  vector  field 
t  =  relative  permittivity  (normally  e  = 

_  F/m) 

E  =  electric  field  intensity  (V/m) 

Cf  =  free  charge  density  (C/mJ) 

H  =  magnetic  field  intensity  (A/m) 
n  =  relative  magnetic  permeability  (nor¬ 
mally  n  =  fu>  =  free-space  magnetic 
permeability  =  4x  x  10"’  H/m) 
w  =  angular  frequency  =  2ir/ 

/  =  frequency  (Hz) 

J  =  current  density  (A/m2) 

j  =  v/=T7 

A  quasi -static  long-wave  length  approximation  is 
employed  in  which  displacement  currents  are  con¬ 
sidered  to  be  negligible  in  comparison  with  conduc¬ 
tion  currents  in  the  model,  and  the  free-space 
wavelength  is  very  long  compared  to  scale  lengths. 


=  free-space  permittivity  =  8.854  xlO'12 


Figure  4.  Transmitting  and  receiving 
coils  over  a  layered  earth. 
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a.  Coaxial  transmitter. 
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b.  Coplanar  transmitter. 

Figure  5.  Primary  electromagnetic  field  direction  for  two 
transmitters.  Arrows  indicate  the  magnitude  dipole  orienta¬ 
tions.  (After  Fraser  1979.) 


This  is  valid  as  long  as  ««  <<  ot,  where  at  is  the  average  conductivity,  and  the  e  of  the  model 
is  on  the  order  of  7  x  10"10  F/m.  Since  typical  AEM  system  bandwidths  are  less  than  50  kHz 
and  scale  lengths  are  normally  much  less  than  the  wavelength  of  50-kHz  radiation,  the  quasi¬ 
static  approximation  is  valid  for  most  AEM  sounding  problems. 

Helicopter  AEM  systems  normally  employ  one  or  more  of  the  three  coil  configurations  de¬ 
picted  in  Figure  1.  There  are  other  coil  arrangements.  A  loop  senses  fields  crossing  its  plane 
as  described  by  Faraday’s  Law  (eq  3)  and  is  only  sensitive  to  field  components  perpendicular 
to  its  plane.  Because  of  this,  the  amplitude  response  of  a  stratified  earth  to  a  vertical  electro¬ 
magnetic  dipole  (whaletail  configuration)  elevated  above  the  surface  is  significantly  different 
from  that  of  a  dipole  in  the  horizontal  plane  (fishtail  or  coaxial  configuration).  The  primary 
magnetic  field  under  the  vertical  coil  is  essentially  horizontal  (Fig.  5).  This  orientation  pro¬ 
vides  excellent  coupling  to  steeply  dipping  bodies.  The  primary  magnetic  field  under  a  hori¬ 
zontal  coil  is  essentially  vertical  and  therefore  couples  well  with  a  horizontal  conductor 
(Fraser  1979).  Figure  6  illustrates  the  coupling  responses  of  vertical  (standard)  coaxial  and 
horizontal  coplanar  (whaletail)  coil  pairs  as  they  move  up  to  and  beyond  vertical  and  inclined 
conductors.  Related  response  variations  obtained  with  the  use  of  a  muiticoil  AEM  system  al¬ 
low  for  detection  and  an  assessment  of  the  depth  and  slope  of  a  conductive  body  or  inter¬ 
face.  Magnetic  dipoles  are  convenient  mathematical  approximations  for  describing  the  re¬ 
sponse  of  loop-type  transmitters,  and  coil  receivers  are  commonly  used  as  sensors  for  the 
magnetic  field.  Maxwell’s  equations  (eq  1-4)  are  now  combined  to  form  the  electromagnetic 
wave  equations  in  E  and  H: 


Figure  6.  Electromagnetic  field  response  difference  between  a  standard  vertical  coaxial 
coil  pair  and  a  whaletail  horizontal  coplanar  coil  pair  moved  across  a  vertical  and  an  in¬ 
clined  dike  conductor.  The  peak  of  the  EM  response  for  the  standard  coil  pair  and  the  trough 
for  the  whaletail  coil  pair  are  shown  to  be  centered  over  the  dike.  (After  Fraser  1979.) 


V1E  =  jog.v>E—  turfE  (5a) 

VW  =  jag.wH- epv’H.  (5b) 

Under  the  quasi-static  approximation  these  expressions  become  diffusion  equations: 

V1E  =  jofiuE  (6a) 

VW  =  joixwH  (6b) 

which  for  a  nonconducting  medium  reduce  to 

vfE^O  (7a) 

V'H  =  0.  (7b) 

Equations  6a  and  b  are  Helmholtz  expressions  with  wave  number  k1  =  -  jnuo. 

Solutions  to  these  equations  can  be  derived  in  terms  of  potentials  suited  to  the  particular 
problem  and  then  differentiated  to  give  the  behavior  of  the  Helds  themselves.  Examples  of 
this  process  can  be  found  in  Wait  (1951).  They  are  often  expressed  as  coupling  ratios  in 
which  the  mutual  coupling  ZT  of  two  coils,  the  magnetic  Tx  dipole  and  the  receiver  coil,  is 
normalized  by  the  “free-space”  coupling  ratio,  for  the  same  coil  configuration,  as  follows 
(Verma  1982): 

Z-r  =  T  (81 


where  ^Rx  =  voltage  at  receiver  coil 
IT%  =  current  in  Tx  coil 
a  =  effective  area  of  coil 
Hr  =  total  magnetic  field  sensed  by  Rx  coil. 

For  the  coplanar  and  coaxial  loops  illustrated  in  Figure  1 ,  the  free-space  coupling  ratio 
values  are,  respectively, 

Z«(cop)  =  +jutu,a,  ch/Aicx'  (9a) 

Zo(coax)  =  -yw/Jofli  a  Jinx' 

=  -2Z0(cop)  (9b) 

where  a,  and  a2  represent  the  effective  area  of  coils  one  and  two,  respectively,  and  x  is  the 
distance  between  coils. 

The  secondary  coupling  ratio  Z,  is  similarly  defined  as 

Z,  =  -jujfia  HJITx  (10) 

where  Ht  is  the  secondary  magnetic  field  at  the  receiver  coil,  which  is  just  the  total  magnetic 
field  minus  the  primary  free-space  field: 


The  secondary  coupling  ratio  just  expresses  the  secondary  magnetic  field  as  a  function  of  the 
primary  field.  It  is  this  ratio,  expressed  in  ppm,  that  was  determined  from  the  AEM  profile 
data. 

Verma  (1982)  reviewed  one  approach  for  computing  these  coupling  ratios  and  provided  a 
convenient  summary  of  the  so-called  Hankel  integrals  T„,  T,  and  Tt  used  in  the  calculation. 
These  integrals  are  defined  as  follows: 

T'(x,t,o,/)=  \  /?,(/, a, X)Xl e-M*. + ».) . J,(\x) dk  (12) 

Ty(x,t_,a,J)  =  [/?,(£,  a.  X)X2e-W'+*->./i(Xx)d\  (13) 

o 

T2(x,l,o,f)  =  [«.(£.  a.  X)Xe-M*-+*.).y,(Xx)dX  (14) 

o 

where  f  =  vector  of  layer  thickness 
a  =  vector  of  layer  conductivity 
X  =  spatial  wave  number 

J»,J,  -  zeroth-  and  first-order  Bessel  function,  respectively 
hit h}  =  height  of  Tx  and  Rx,  respectively  (Fig.  4) 

R,  =  reflection  coefficient,  defined  recursively  by  Koefoed  (1972)  as 


where  Ra  a  =  0 

Vx  .  sfr’  +  kf 
*l  =  sl-junOi 
*U  =  (K,-Kk)/(KI+Kk). 

The  horizontal  coplanar,  vertical  coplanar  and  vertical  coaxial  loop  secondary  coupling 
ratios  are  now,  respectively 

£=*T.  (16) 

i  =  (17) 

J^=j t'cr.-xr,)^.  (is) 

The  technique  reviewed  by  Verma  (1982)  and  first  used  by  Ghosh  (1971)  is  commonly  called 
the  "linear  digital  filter”  method.  It  permits  the  computation  of  rapid,  stable  estimates  of 
the  electromagnetic  response  of  a  layered  earth  model.  The  technique  works  very  well  and 
serves  as  an  excellent  basis  for  a  general  forward  and  inverse  modeling  technique.  However, 
for  our  AEM  interpretation  requirements,  his  approach  represents  “overkill.”  A  much 
faster  method  for  performing  the  integrals  (eq  12-14)  was  used,  based  on  Gauss-Laguerre  in¬ 
tegration  of  the  Hankel  transforms.  Gauss-Laguerre  is  a  standard  technique  of  numerical 
analysis.  The  integrals  are  performed  five  to  ten  times  faster  than  for  the  shortest  acceptable 
digital  filter  evaluated.  Since  the  inversion  process  spends  most  of  its  time  on  the  evaluation 
of  the  Hankel  transform  integrals,  the  inversion  is  also  speeded  up  by  a  factor  of  five  to  ten. 

The  inversion  method  we  used  to  analyze  the  AEM  data  is  based  on  a  standard  numerical 
technique  of  least-squares  nonlinear  regression  with  singular  value  truncation  and  damping, 
optimized  for  ice  and  water  depth  interpretation.  All  four  frequencies  (eight  channels)  of 
AEM  data  are  used  together  in  the  inversion  process.  It  is  a  generalized  and  enhanced  version 
of  the  Gauss-Newton  and  Marquardt  methods  (Marquardt  1963).  Consider  the  computed  re¬ 
sponse  C7  of  a  layered  earth  model,  a  known  nonlinear  function  of  the  model  parameter  a,  £ 
and  p: 

C  =  F(a,t_,p)  (19) 

where  p  represents  the  system  parameters  such  as  coil  separation,  frequencies,  altitude, 
pitch,  roll,  etc.  However,  Fean  be  approximately  linearized  by  a  Taylor  series  expansion  in 
the  model  parameters  a  and  £.  Thus,  if  X_  represents  the  desired  model  parameters  (including 
perhaps  both  a  and  £),  F  can  be  expressed  as  the  infinite  series 

dF(X »)  d2F(Xt) 

F(X-Xt)  =  F(X,)+  ~~  •  (X~X,)+  -j=-  .  (X-X,y  +  .  .  .  (20) 

where  Xj>  represents  a  set  of  parameters  for  the  “starting  model”  in  the  inversion.  With  the 
starting  model  X_ ,  “close”  to  the  true  model  Xt,  the  series  can  be  approximated  by  trunca¬ 
tion  to  the  first-order  term 


dF(Xo) 

F(X-X,)  s  F(A>)+  •  (X-Xt). 


(21) 


7 


This  “linearization”  process  is  now  set  up  in  terms  of  the  residual  response  AC  =  0  -  and 
a  parameter  correction  vector  AX  =  X_- X*,  where  X»  is  now  the  base  model  for  this  step 
of  the  inversion  process,  0  is  a  set  of  observations  and  C_  =  F(Xt).  Then,  from  eq  21  we  ob¬ 
tain 


dF(X,) 

AC  =  jjf  •  Ap  (22) 

in  which  dF/dX_  is  the  Jacobian  matrix  =  df{/dX},  where  f^X)  =  Q.  This  matrix  is  nor¬ 
mally  non-negative  definite  and  cannot  always  be  inverted,  even  with  infinite-precision  arith¬ 
metic,  because  of  the  zero  eigenvalues  that  it  may  possess.  To  get  around  this  problem,  ap¬ 
proximate  inverses  are  used  (Jupp  and  Vozoff  1957,  Holladay  1980).  Briefly  the  matrix  is  de¬ 
composed  into  two  eigenvector  matrices  and  a  set  of  singular  values.  This  process  is  called 
the  singular  value  decomposition  by  Lancsoz  (1961)  and  is  expressed 

J  =  U  A  V'  (23) 

where  J  =  Jacobian  matrix 
U  =  data  eigenvector 
V  =  parameter  vector 

A  =  diagonal  matrix  of  singuar  values  for  J. 

In  this  process,  rows  and  columns  of  U  and  V  corresponding  to  zero  and  near-zero  singular 
values  are  discarded,  and  damping  is  applied  to  the  remaining  singular  values  as  required. 
The  damped  matrix  is  then  inverted  to  give  a  solution 

AX  =  VA-'U'AC  (24) 

which  is  used  to  correct  the  parameter  vector  X.  Then  a  new  forward  modelFf^O  is  com¬ 
puted,  a  new  residual  AC  is  formed,  and  the  process  is  repeated  until 

IacI  =  Vac'*  ac  (25) 


falls  within  acceptable  limits. 

The  starting  model  is  fairly  important  to  the  efficiency  of  the  process.  Indeed,  if  the  start¬ 
ing  model  is  too  “far”  from  the  true  solution,  the  process  may  not  converge.  To  avoid  this 
problem,  a  table  was  used  to  prepare  a  good  starting  model. 

The  above  analysis  will  always  allow  determination  of  the  distance  from  the  AEM  bird  to 
the  seawater  surface,  provided  the  bird  is  flown  within  a  prescribed  elevation  window.  This 
elevation  window  is  dictated  by  the  AEM  system  response,  sensitivity  and  dynamic  range, 
and  is  generally  20-50  m  above  the  surface.  Flying  the  bird  below  about  20  m  may  cause  the 
electromagnetic  Held  strength  from  seawater  to  saturate  the  receiver.  This  can  be  overcome 
by  reducing  the  receiver  sensitivity  but  currently  not  in  real  time.  Flying  the  bird  higher  than 
about  50  m  can  cause  the  electromagnetic  field  from  the  seawater  to  be  too  weak  to  measure 
accurately.  As  an  examination  of  the  EM  response  equations  will  show,  the  electromagnetic 
bathymetry  problem  will  not  always  allow  for  accurate  inversion.  There  is  the  problem  of 
screening  by  the  conductive  seawater  and,  to  a  much  lesser  extent,  by  conductive  sea  ice.  As 
the  electromagnetic  wave  penetrates  the  ice  and  water,  it  is  attenuated  and  phase-rotated.  Re¬ 
flections  from  the  water  surface  and  from  the  seabed  are  likewise  attenuated  progressively  by 
greater  ice  and  water  thickness  and  by  flight  height,  should  it  increase.  Practical  limitations 
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Figure  7.  Nomogram  for  obtaining  repre¬ 
sentative  electromagnetic  skin  depth  vs  fre¬ 
quency  and  conductivity.  (After  Won  1980.) 


on  the  frequencies  used  in  the  AEM  bird  therefore  limit  the  maximum  depth  of  water  that 
can  be  interpreted  at  a  given  error  level,  as  well  as  the  maximum  water  depth  at  which  the 
conductivity  of  the  seafloor  can  be  estimated.  For  practical  purposes,  and  with  error  levels 
associated  with  the  AEM  system  used,  the  error  is  less  than  5%  of  the  water  depth  down  to 
one  skin  depth.  The  probable  error  in  water  depth  for  the  AEM  system  used  is  ±0.3,  ±  1 .0, 
and  +2  and  -2  m  for  water  depths  of  10,  13  and  20  m,  respectively.  As  a  general  estimate 
the  error  is  about  10%  of  the  seawater  depth  at  about  1.3  skin  depths  for  the  system  used. 

The  electromagnetic  skin  depth,  or  the  depth  of  penetration,  is  mainly  determined  by  the 
AEM  system  frequency  and  the  ground  (seawater)  conductivity.  For  most  situations  the  skin 
depth  D  defines  the  limiting  depth  at  which  reasonable  AEM  sounding  estimates  can  be 
made.  The  electromagnetic  skin  depth  equals  -fl/wyo  =  303  \JWof.  For  seawater  the  sound¬ 
ing  depth  can  be  about  1 .3  skin  depths.  The  nomogram  in  Figure  7  gives  the  relationship  be¬ 
tween  AEM  system  source  frequency,  ground  conductivity  and  skin  depth  for  a  wide  range 
of  materials  and  sounding  frequencies.  For  seawater  the  skin  depth  is  better  illustrated  in 
Figure  8.  Since  Beaufort  Sea  water  during  the  winter  has  a  conductivity  of  about  2.3  S/m, 
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Figure  8.  Electromagnetic  skin  depth  in  saline  water  vs  anten¬ 
na  frequency  and  water  conductivity. 


IP,  In  Phose  (ppm) 


Figure  9.  Example  of  an  Argand  diagram  used  for  determining 
the  AEM  bird  elevation  above  (and  the  apparent  conductivity  of) 
a  conductive  half-space.  The  diagram  models  the  response  of  a  hori¬ 
zontal  coptanar  coil  pair,  separated  about  6.5  m. 


the  effective  sounding  depth  of  our  AEM  system  at  527.5  Hz  is  between  about  13.5  (1  skin 
depth)  and  21  m  (1 .5  skin  depth).  The  conductivity  of  the  surface  layer  of  Beaufort  Sea  water 
in  summer  will  generally  be  lower  because  of  melting  of  the  sea  ice  canopy.  Under  this  condi¬ 
tion  the  AEM  system  sounding  depth  will  be  proportionately  deeper. 

With  the  aid  of  an  Argand  diagram  (Fig.  9)  formulated  on  the  basis  of  a  given  Tx-Rx  coil 
spacing  and  orientation,  the  height  of  the  bird  above  the  condutive  surface  and  the  apparent 
conductivity  of  the  conductive  layer  can  be  estimated.  For  example,  if  the  in-phase  and 


y> 


y>y>y>y>y»ytytyr, 


quadrature  amplitudes  at  the  receiver  coil  are  2000  and  850  ppm,  respectively,  then  the  inter¬ 
cept  of  the  related  lines  in  Figure  9  indicates  that  the  bird  is  about  26  m  above  the  conductor, 
or  in  our  case  the  seawater,  and  that  the  response  parameter  is  5000.  It  follows  that  if  the  Tx 
coil  was  operating  at  a  frequency  of  2000  Hz,  the  apparent  seawater  conductivity  would  be 
5000/2000  =  2.5  S/m.  Through  appropriate  multilayer  analyses  as  outlined  above,  the 
depth  to  and  conductivity  of  multiple  conductive  layers  can  be  estimated. 

Normally  for  AEM  interpretation,  three  “layers”  are  considered  in  the  inversion  process: 
ice,  water  and  bottom.  The  data  collected  in  1985  were  not  adequate  to  enable  determination 
of  ice  and  bottom  sediment  conductivity,  so  the  parameters  determined  by  inversion  were  ice 
thickness,  water  depth  and  water  conductivity.  Ice  and  sea  bottom  conductivities  were  em¬ 
pirically  fixed  at  reasonable  values  based  on  preliminary  data  interpretations  and  on  compar¬ 
isons  of  these  findings  to  available  ground  truth  at  specific  sites.  The  conductivities  for  both 
the  sea  ice  and  the  sediment,  typically  about  0.02  S/m,  are  very  low  compared  to  the  sea¬ 
water. 

A  reason  for  reducing  the  number  of  parameters  interpreted  is  to  increase  the  volume  of 
EM  data  “available”  for  estimating  the  snow  and  ice  thickness.  For  some  interpretations  the 
seawater  conductivity  was  fixed  at  2.5  or  2.6  S/m  to  observe  the  effect  on  the  snow  and  ice 
thickness  and  water  depth  determinations.  Of  the  examples  given  in  this  paper,  only  for  line 
6L7  was  the  seawater  conductivity  fixed.  This  helped  reduce  noise  or  scatter  in  the  interpret¬ 
ed  snow  and  ice  thicknesses  and  water  depths.  Fixing  the  seawater  conductivity  at  a  reasona¬ 
ble  value  does  not  affect  the  snow  and  ice  thickness  determination  very  much  since  the  highly 
conductive  water  yields  a  strong  AEM  response,  which  allows  a  good  determination  of  snow 
and  ice  thickness  even  if  the  water  conductivity  used  is  not  very  accurate. 

If  system  noise  levels  are  1  ppm,  water  conductivity  and  distance  to  a  horizontal  water  sur¬ 
face  should  be  estimated  to  about  17s.  Sea  ice  thickness  should  be  estimated  to  within  a  few 
percent  for  thick  plate  or  level  ice.  However,  because  the  laser  profiler  used  had  an  accuracy 
on  the  order  of  ±  10  cm,  ice  thickness  estimates  will  vary  accordingly. 

Ice  conductivity  may  also  be  estimated  but  only  with  accuracy  if  the  high-frequency  coils 
operate  well  above  the  32-kHz  frequency  available  for  this  study.  A  50-kHz  frequency  will  be 
used  in  our  1987  field  program.  However,  Liu  and  Becker  (1987)  indicated  that  an  order  of 
magnitude  higher  frequency  may  be  needed  to  determine  the  conductivity  of  sea  ice. 


CALIBRATION  OF  AIRBORNE  ELECTROMAGNETIC  DATA 

Typical  AEM  systems  measure  the  amplitude  of  the  secondary  field  at  the  receiver  as  a 
fraction  of  the  primary  field  (from  the  transmitter)  in  parts  per  million  (ppm).  The  phase  of 
the  measured  signal  can  also  be  shifted  relative  to  the  transmitted  field.  The  amount  of  this 
phase  shift  contains  as  much  information  as  the  overall  amplitude  of  the  signal.  However, 
instead  of  measuring  the  amplitude  and  phase  of  the  signal,  the  response  is  usually  measured 
as  the  amplitude  of  each  of  the  in-phase  and  out-of-phase  (quadrature)  signals. 

Calibration  and  phasing  of  an  AEM  system  consist  of  determining  the  size  of  these  in- 
phase  and  quadrature  responses  from  the  actual  voltages  measured  and  recorded  by  the  data 
acquisition  system.  Calibration  thus  yields  a  factor  in  parts  per  million  (of  the  primary  field 
at  the  receiver)  per  volt  (ppm/V),  which  is  later  multiplied  by  incoming  survey  data  (in  volts) 
to  obtain  the  correct  response  in  ppm. 

For  the  Prudhoe  Bay  survey,  phasing  and  external  calibrations  were  carried  out  on  the 
ground  using  techniques  which  are  standard  for  mineral  exploration  surveys.  Phasing  en¬ 
sures  that  the  signal  used  as  a  time  reference  in  measuring  phase  shift  has  the  correct  phase  it¬ 
self.  Once  calibration  and  phasing  have  been  completed,  factors  for  both  the  in-phase  and 


quadrature  data  are  determined  by  applying  a  known  secondary  field  to  the  receiver  and 
measuring  the  system  response  in  volts. 

However,  this  standard  calibration  technique  for  AEM  data  is  not  sufficiently  accurate  for 
the  more  stringent  requirements  of  ice  measurement  or  bathymetry.  Therefore,  calibrations 
had  to  be  performed  on  a  more  empirical  basis  for  this  theory.  Initially,  forward  models 
were  computed  to  obtain  theoretical  in-phase  and  quadrature  data  at  sites  where  the  ice 
thickness  and  water  depth  and  conductivity  were  "known.”  The  model  data  (in  ppm)  were 
then  compared  to  the  real  EM  data  (in  ppm)  measured  at  each  ground  truth  site.  Changes  to 
the  initial  calibration  factors  were  then  made  to  yield  an  improved  calibration.  This  tech¬ 
nique  is  directly  comparable  to  the  use  of  a  "bar  check”  for  echo  sounder  calibrations.  How¬ 
ever,  the  use  of  ground  truth  for  calibration  restricts  the  accuracy  of  the  AEM  system  results 
to  the  accuracy  of  the  ground  truth  data.  It  also  requires  that  ground  truth  be  obtained  on 
every  survey,  e.g.  flying  down  a  large  unfrozen  lead  located  in  deep  water.  When  a  good  cali¬ 
bration  is  obtained  for  an  AEM  system  in  the  field,  it  has  been  found  to  remain  very  consis¬ 
tent  from  day  to  day. 

Zero  levels  needed  in  the  analysis  of  baseline  drift  removal  (a  practice  common  to  a  wide 
variety  of  high-resolution  measurement  techniques)  are  obtained  by  flying  the  system  to  an 
altitude  of  up  to  500  m.  At  this  elevation,  there  is  no  measurable  response  from  the  seawater. 
Zero  level  measurements  were  made  before  and  after  each  flight  line.  Drift  is  normally  as¬ 
sumed  to  be  linear  with  time  between  these  zeroing  measurements,  although  in  practice  some 
nonlinearity  is  sometimes  evident.  It  is  thought  that  this  drift  is  caused  by  small  changes  in 
the  temperature  of  certain  electronic  components  during  a  survey  flight. 

FIELD  SURVEY 

The  AEM  system  flights  were  made  over  a  variety  of  ice  formations.  These  include  first- 
year  sea  ice  of  several  thicknesses  and  first-  and  second-year  pressure  ridges. 

Newly  formed  lead 

To  assess  system  response  and  to  aid  in  system  calibration,  a  flight  was  made  down  the 
lead  shown  in  Figure  10.  The  lead  was  in  the  process  of  freezing  over.  The  near  end  of  the 


Figure  10.  Skim  ice  on  a  newly  formed  lead  profiled  with  the  AEM  system. 


lead  had  about  1  cm  of  ice  cover,  while  the  far  end  was  still  open  water.  For  our  purposes, 
the  lead  represented  a  smooth,  conductive  water  surface.  Indeed,  the  AEM  system  responses 
indicated  this  by  profiling  a  level  surface  over  this  feature.  An  average  water  conductivity  of 
2.8  S/m,  about  l(Wt  too  high,  was  also  determined  from  the  AEM  data. 
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Refrozen  lead 

A  profile  made  over  a  uniformly  thick  refrozen  lead  with  a  3-cm  snow  cover  (Fig.  1 1)  is 
shown  in  Figure  12,  and  the  data  related  to  the  profile  are  given  in  Table  A1 .  The  calibrated 
and  drift-corrected  AEM  data  for  this  profile  are  given  in  Figure  13.  The  profile  in  Figure  12 
indicates  a  non-uniform  lead  ice  thickness.  This  variation  is  not  due  to  snow  cover  variations 
but  is  believed  to  be  due  to  system  noise  and  drift  problems,  to  the  poor  accuracy  of  the  laser 
profilometer,  as  previously  mentioned,  and  to  bird  pitch  and  roll  variations.  The  latter  could 
not  be  fully  accounted  for  in  the  bird  pendulum  data.  Nevertheless,  the  average  ice  thickness 


Figure  11.  Refrozen  lead  profiled  with  the  AEM  system. 


Figure  12.  Snow  and  ice  thickness  and  under-ice  water  depth  for  a  flight 
over  a  refrozen  lead  located  between  about  fid.  no.  2880  and  3000. 
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Figure  13.  AEM,  laser  altimeter  and  altitude  sensor  data  for  flight  over  refrozen  lead  (Fig.  11  and  12). 

The  graphs  give  the  bird  roll  and  pitch  in  degrees,  the  laser  altitude  in  meters,  and  the  16,020-,  4160- ,  930-  and 

530-Hz  frequency  coil  in-phase  (IV)  and  quadrature  phase  (Q)  response  in  ppm. 

for  the  profile,  given  in  Table  Al,  is  0.80  m.  This  is  remarkably  close  to  the  measured  ice 
thickness  of  0.75  m.  This  result  is  extremely  encouraging  and  was  achieved  after  significant 
improvements  were  made  to  the  computer  program  to  account  for  system  noise  and  bird 
pitch  and  roll.  The  first  analysis  of  the  data  gave  a  lower  average  ice  thickness  of  0.65  m 
(Kovacs  et  al.  1987).  An  improved  vertical  accelerometer  and  pitch  and  roll  sensor,  a  reduc¬ 
tion  of  system  noise  and  a  more  accurate  laser  profilometer  should  further  improve  the  pro¬ 
file  results.  The  estimated  seawater  conductivity  ctw,  listed  in  Table  Al,  varied  from  2.35  to 
2.49  S/m.  The  average  was  2.45  S/m,  which  is  in  very  good  agreement  with  measured  values 
of  2.5  S/m.  Again,  the  first  analysis  of  the  data  gave  an  average  ow  of  3.0  S/m,  or  about  0.5 
S/m  higher  than  the  measured  values.  The  seabed  was  also  profiled  at  this  site.  This  was  not 
anticipated  and  therefore  no  direct  sounding  measurements  were  made  to  verify  the  AEM 
depth  determinations.  The  variation  in  the  bottom  profile  is  not  real,  however,  and  is  prob¬ 
ably  due  to  the  factors  discussed  above. 

Second-year  floe 

A  second-year  ice  floe  was  selected  for  AEM  profiling.  The  floe  had  a  ridge  extending 
from  one  side  to  the  other.  Survey  lines  were  established  on  this  floe  by  an  Exxon  research 
team.  The  lines  were  spaced  23  m  apart,  ran  perpendicular  to  the  ridge  and  were  250  m  in 
length.  We  used  two  of  these  lines  and  ran  a  third  line  between  them.  This  resulted  in  a  test 
area  250  m  long  with  three  parallel  lines  about  11.5  m  apart  (Fig.  14).  Along  each  line  the 
snow  and  ice  thickness  was  measured  by  drilling  holes  at  about  7.5-m  intervals  except  at  the 


e 


i'  jV  .  ■ 


•  *  A 


&2*  EWAdSSS 


V.' 


Jjmt ;  v  :.;^  „, ■*&?•''  **  *  '.&- 

W«  '•  '  <  -J^r* 4  *L  ' 

• :  '  ' -K  ^  4  **  -  •  ’ 

I*~  .  y-^w  &*■'&*'  "  &*••"••  '  ^  v 

vw'  *  'fc  ' 

'*w  TOv  -»  “•  ■ 

ste-  "  jp$/  _ 

r^v 

'  *3^ 

B*  .*  ■ 


,»;  .y -v 
;  .P.-*rP: 


7L*>, 


f«’ 


-.?**■  v 


■«,.•,  .**■ '  *  ,r 


.  * , n 

>  jf-  ’IT 

jd  4>"-^  &■ 

'f  •  ~  *.V 

*  T*  , 


■  *>*—■  '  i  > 

Figure  IS.  Aerial  view  of  second-year  sea  ice  flow  and  location  of  survey  lines  A,  B 
and  C. 


ridge,  where  the  spacing  was  about  3.75  m  on  lines  A  and  B.  At  30-m  intervals  along  each 
line,  we  placed  markers  that  could  be  seen  from  the  air.  Figure  15  is  an  aerial  view,  looking 
east,  of  the  undulating  floe  relief,  the  ridge  and  the  relative  location  of  the  three  parallel 
lines.  Two  men  can  be  seen  standing  on  line  A  (near  side  of  the  ridge)  and  the  person  on  the 
trail,  arrow  at  upper  right,  gives  a  perspective  of  scale.  The  large  dark  object  on  the  opposite 
side  of  the  ridge  is  a  1.2-m-square,  2.2-m-high  plastic  building. 

The  snow  and  ice  station  thickness  measured  along  each  survey  line  is  given  in  Table  1  and 
is  shown  graphically  in  Figure  16  along  with  the  average  snow  and  ice  thickness  for  all  three 
lines.  The  cross  section  of  the  ridge  between  stations  13  and  19  for  each  survey  line  is  given  in 
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Figure  16.  Snow  and  ice  thickness  along  survey  lines  A,  B  and  C,  and 
average  snow  and  ice  thickness  for  the  three  lines. 

Figure  17.  The  cross  sections  show  that  the  ridge  sail  is  just  under  2  m  high,  the  deepest  part 
of  the  keel  is  offset  from  the  center  of  the  sail,  the  keel  depth  and  geometry  are  different  for 
each  line,  and  the  snow  cover  is  thickest  in  the  area  of  the  ridge,  which  acts  as  a  “snow 
fence.”  The  snow  and  ice  thickness  profiles  (Fig.  16)  show,  as  expected,  that  variations  in 
thickness  exist  between  the  profiles  at  the  same  station  locations.  Again,  the  thickest  snow 
and  ice  was  located  at  the  ridge,  which  is  situated  between  stations  IS  and  17.  The  average 
snow  and  ice  thickness  for  profiles  A,  B  and  C  was  3.56,  3.58  and  3.74  m,  respectively,  for 
an  overall  average  of  3.62  m. 

The  original  intent  was  to  fly  the  AEM  bird  down  each  survey  line  and  note  each  ice  sur¬ 
face  station  location  on  the  data  file  as  the  bird  passed  over  the  marker.  Then  a  comparison 
of  the  AEM  ice  thickness  data  to  the  drill  hole  data  was  to  be  made.  This  was  not  possible  for 
reasons  given  below.  Flights  were  made  from  west  to  east,  nearly  into  the  prevailing  wind. 
The  first  flight  pass  quickly  revealed  that  it  would  be  extremely  difficult  to  fly  the  bird,  sus¬ 
pended  30  m  below  the  helicopter,  down  such  an  exact  linear  course.  Given  this  difficulty,  an 
effort  was  made  to  fly  the  bird  as  close  as  possible  down  the  center  line  of  the  three-line 
track. 

The  snow  and  ice  thickness  as  obtained  from  two  passes,  flights  6L3  and  6L4,  with  the 
16-kHz  coil  in  the  bird  as  shown  in  Figure  18.  The  thickness  for  three  passes  with  the  32-kHz 
coil,  flights  8L2, 8L3  and  8L4,  is  shown  in  Figure  19.  The  data  associated  with  each  flight  are 
listed  in  Tables  A3-A7. 
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Flight  6L4  was  made  at  the  lowest  flight  speed,  during 
which  51  thickness  measurements  were  made  down  the 
250-m  track.  This  was  followed  by  flight  8L2,  during 
which  34  thickness  measurements  were  made.  Flights 
6L3,  8L3  and  8L4  were  made  at  about  the  same  speed, 
during  which  26,  26  and  28  thickness  measurements,  re¬ 
spectively,  were  made. 

Snow  and  ice  thickness  profiles  from  flights  6L4  and 
8L2  (Fig.  18b  and  19)  show  a  waviness  not  at  all  charac¬ 
teristic  of  the  measured  thickness.  During  these  two 
slow  runs  the  bird  was  not  stable.  It  moved  from  one 
side  of  the  flight  path  to  the  other.  The  nose  of  the  bird 
would  turn  slowly  left  and  then  move  beneath  the  heli¬ 
copter  to  the  left  side;  then  the  nose  would  swing  right, 
and  the  bird  would  slide  to  the  right  like  a  pendulum.  This  movement  could  not  be  fully  ac¬ 
counted  for  in  the  analysis  of  the  data  and  is  believed  to  be  the  cause  of  the  undulations  in 
the  snow  and  ice  thickness  profile.  The  profiles  obtained  at  the  higher  speed  are  more  repre¬ 
sentative  of  the  average  measured  snow  and  ice  thickness  profile  in  Figure  16  in  that  there  is  a 
gradual  thickening  of  the  ice  in  the  area  of  the  ridge.  Nevertheless,  the  AEM  profiles  do  not 
clearly  reveal  the  ridge  nor  do  they  show  the  maximum  ice  thickness  measured  by  drilling. 

The  AEM  data  do  not  show  the  thick  ridge  ice  because  of  the  large  footprint  and  therefore 
the  surface  area  over  which  the  water  surface  was  integrated  into  each  AEM  distance  deter¬ 
mination.  Remember,  the  AEM  system  determines  the  distance  from  the  bird  to  the  water 
surface,  which  has  depressions  in  it  caused  by  ice  bottom  relief  variations.  Our  preliminary 
assessment  is  that  the  AEM  system  footprint  diameter  is  about  1.25  times  the  bird’s  height. 
This  was  determined  by  averaging  the  drill-hole-measured  snow  and  ice  thickness  data  over  a 
window  until  the  average  thickness  agreed  with  the  AEM  ice  thickness  interpretation  at  the 
center  of  the  window.  This  footprint  size  has  recently  been  verified  by  A.  Becker  and  G. 
Liu,*  who  made  an  analysis  similar  to  ours  and  obtained  comparable  results.  Liu  and  Becker 
also  made  a  worst-case  assessment  based  on  the  expected  current  density  on  a  flat,  perfectly 
conducting,  seawater  interface,  and  then  determined  the  relative  area  that  would  account  for 
90%  of  the  observed  secondary  magnetic  field.  Their  conclusion  was  that  the  footprint  of  an 
AEM  system  is  about  1.4  times  the  bird  height  for  a  horizontal  double-dipole  coil  system 
(with  vertical  coils).  Therefore,  the  AEM  distance  to  the  ice/water  interface  is  averaged  over 
a  relatively  large  area  of  undulating  sea  surface  relief.  This  effect  smooths  out  the  snow  and 
ice  thickness  variations  in  the  drill-hole-measured  thickness  profile  (Fig.  16). 

Smoothing  the  snow  and  ice  relief  by  the  AEM  system  would  clearly  be  of  concern  to  those 
interested  in  detecting  and  measuring  pressure  ridge  cross  sections  and  depths.  Nevertheless, 
this  does  not  weaken  the  case  for  AEM  measurement  of  snow  and  ice  thickness.  According 
to  Becker  and  Liu,  “it  only  strengthens  the  case  for  developing  better  methods  of  data  inter¬ 
pretation.”  However,  for  measuring  the  mean  snow  and  ice  thickness  of  the  arctic  pack  ice, 
AEM  profiling  should  prove  very  useful.  This  may  be  examined  in  Table  1,  which  lists  the 
average  snow  and  ice  thickness  for  each  of  the  five  flights  down  the  track  on  the  second-year 
ice  floe,  as  well  as  the  average  seawater  conductivity.  The  drill-hole  measurements  gave  an 
average  snow  and  ice  thickness  of  3.62  m  for  the  three  survey  lines.  The  highest  average 
AEM  snow  and  ice  thickness,  4.07  m,  is  for  flight  8L2  and  the  lowest,  3.07  m,  is  for  flight 
8L3.  These  values  are  12.4%  higher  and  15.2%  lower  than  the  average  directly  measured 


Table  1.  Average  snow  and  ice 
thickness  f,  and  seawater  con¬ 
ductivity  a,  for  each  flight  over 
the  second-year  sea  ice  floe. 


Flight 

no. 

Avg.  la 
(m) 

Avg.  aw 
(S/m) 

8L2 

4.07 

2.59 

8L3 

3.07 

2.35 

8L4 

3.55 

2.46 

6L3 

3.20 

2.48 

6L4 

3.45 

2.47 

Avg.  all 

3.47 

2.47 

•  Personal  communication.  University  of  California  at  Berkeley. 


value.  From  this  analysis,  the  AEM  system  provided  a  representative  mean  snow  and  ice 
thickness  to  within  about  15%  of  that  measured  by  direct  drilling. 

Other  factors  may  have  affected  the  results.  The  difficulty  of  flying  the  AEM  bird  down 
the  center  of  the  track  for  which  representative  snow  and  ice  thickness  data  existed  was  cer¬ 
tainly  one.  Another  problem  is  the  three-dimensional  form  of  the  ice  features.  The  106  drill¬ 
hole  measurements  made  along  the  track  may  not  have  been  sufficient  to  fully  define  the 
mean  ice  thickness,  and  the  algorithms  used  to  interpret  the  AEM  data  were  not  refined 
enough  to  handle  3-D  ice  features  (water  depressions).  The  inversion  technique  is  being  im¬ 
proved,  and  the  bird  coil  arrangement  will  be  changed  to  allow  better  detectability  and  defi¬ 
nition  of  the  ridge  keel/seawater  interface. 

The  average  seawater  conductivities  listed  in  Table  1  were  in  good  agreement  with  the 
measured  value  of  2.5  S/m. 


Grounded  floeberg 

In  the  early  winter  of  1983  a  severe  storm,  with  winds  out  of  the  southwest,  hit  the  Beau¬ 
fort  Sea  coast  in  the  area  of  Prudhoe  Bay,  Alaska.  The  high  winds  broke  up  the  sea  ice  over 
much  of  the  continental  shelf  and  drove  this  ice  and  the  offshore  pack  ice  westward.  During 
this  event,  four  deep-keeled  “multi-year”  ice  floes  became  grounded  on  Cat  Shoals.  Sea  ice 
floes  driven  against  these  anchor  points  failed.  At  first  the  ice  piled  on  top  of  the  grounded 
ice  floes,  adding  weight  and  increasing  their  sliding  resistance.  The  ice  continued  to  break  for 
some  time.  The  ice  rubble  increased  in  height  and  then  began  to  grow  upwind.  When  the 
storm  had  driven  the  last  of  the 


first-year  ice  out  to  sea,  four  very 
large  grounded  ice  formations  re¬ 
mained.  These  features  were  about 
0.5  to  0.75  km  long  and  up  to  100  m 
wide.  The  location  of  these  forma¬ 
tions  on  Cat  Shoals,  as  determined 
by  global  navigation,  is  given  in  Fig¬ 
ure  20.  The  bathymetric  high  points 
on  Cat  Shoals  are  typically  less  than 
10  m  below  the  surface,  but  shoaling 
to  8  m  at  one  location  is  shown  on 
one  bathymetric  chart  of  the  area. 
Ice  formation  D  seems  to  be  ground¬ 
ed  in  waters  deeper  than  12  m,  while 
formations  A-C  appear  to  have  ini¬ 
tiated  at  a  grounding  point  in  water 
12  m  or  less  in  depth.  The  water 
depths  at  the  southwest  ends  of  for¬ 
mations  A-D  were  measured  by  lead 
line  and  found  to  be  1 1 .5, 1 1 .6, 12.5 
and  12.5  m,  respectively.  The  two 
largest  formations  were  B  and  C. 
These  features  each  had  an  ice  ridge 
over  15  m  high  (Fig.  21).  Higher 
grounded  sea  ice  rubble  formations 
do  occur.  One  that  formed  on  a 
shoal  about  23  km  north  of  Oliktok 
Point  was  22  m  high  (Fig.  22). 

The  four  grounded  ice  formations 
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on  Cat  Shoals  did  not  break  up  or  Figure  20.  Location  of  grounded  ice  formations  on 


lose  much  ice  by  calving  during  the  Cat  Shoals,  and  of  Loon  Shoal,  Sandpiper  and  North 


1984  summer.  Figure  23  shows  the  Star  islands,  and  AEM  sounding  track  6L7. 


Figure  21.  View  of  south  side  of  grounded  rubble  formation  C  and  its  highest  ridge. 


Figure  22.  Extremely  high  (22  m)  rubble  formation  found  grounded  in  1984 
on  a  shoal  north  of  Oliktok  Pt.,  Harrison  Bay,  Alaska,  at  about  70°42. 7  'N, 
149°44'W. 

configuration  of  formations  B,  C  and  D  in  late  August  1984,  and  Figure  24  shows  the  con¬ 
figuration  of  formations  B  and  C  on  27  March  1985.  Having  survived  the  summer  melt 
season,  these  formations  may  now  be  referred  to  as  grounded  floebergs.  At  this  time 
floeberg  B  was  about  650  m  long  and  80  m  wide,  and  floeberg  C  is  some  600  m  long  and  up 
to  90  m  wide. 

Stereo  photography  techniques  were  used  to  determine  the  topography  of  floeberg  C  along 
the  parallel  lines  shown  in  Figure  25.  The  lines  are  15  m  apart.  The  highest  point  on  the  floe- 
berg  was  about  13.5  m  high  on  line  4  at  the  location  of  the  circle.  Figure  26  is  a  low-level 
aerial  view  of  this  high  point,  and  Figure  27  is  a  view  from  the  ice. 


C  and  D  in  August  1984. 


Figure  25.  Floeberg  C,  showing  the  location  of  the  highest  ridge  and  the  lines  along 
which  the  topography  was  determined  via  stereo  analysis. 


Figure  27.  North  of  floeberg  C  and  its  highest  ridge.  Note  the  vertical  face  where  ice 
had  calved  off  during  the  previous  summer;  this  ice  face  is  voidless. 
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Figure  28.  AEM  profile  results  from  the  flight  over  floeberg  C.  Ridges  A-E  and  num¬ 
bered  subpeaks  correspond  to  similar  features  in  Figure  29. 


An  AEM  sounding  run  was  made  down  each  of  the  floebergs.  The  snow  and  ice  thickness 
and  under-ice  water  depth  profiles  for  floeberg  A  were  presented  in  Kovacs  et  al.  (1987). 
Here  we  give  AEM  sounding  results  for  floeberg  C.  The  flight  down  floeberg  C  was  made  by 
lining  up  with  the  high  point  on  the  floeberg.  This  path  would  be  nearly  that  along  line  4. 
The  AEM  snow  and  ice  thickness  and  under-ice  water  depth  profiles  are  shown  in  Figure  28, 
and  the  tabulated  data  are  given  in  Table  A8  along  with  bird  pitch,  roll  and  flight  elevation 
and  seawater  conductivity.  The  AEM  results  indicate  there  is  an  average  of  4.8  m  of  water 
under  the  floeberg.  This  is  not  possible  since  the  floeberg  is  solidly  grounded  on  the  shoal. 
The  cause  of  this  ambiguity  is  believed  to  be  the  existence  of  seawater  within  the  submerged 
ice  keel  block  structure.  This  is  not  unreasonable.  The  AEM  system,  because  of  the  large 
footprint,  may  also  have  sensed  the  seawater  off  to  the  side  of  the  floeberg.  In  addition,  be¬ 
cause  the  floeberg  is  three-dimensional,  the  layered  model  used  in  the  data  analysis  may  not 
have  been  adequate  to  define  this  imperfect  transition  between  the  ice  block,  the  water  and 
the  seabed.  Further  refinement  of  the  AEM  data  interpretation  technique  is  progressing.  A 
better  solution  to  this  apparent  definition  problem  is  anticipated. 
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Figure  29.  Snow  and  ice  elevation  profiles  as  determined  by  stereo  analysis  techniques  along  lines  3,  4 
and  5  in  Figure  25.  Ridges  A-E  and  numbered  subpeaks  correspond  to  similar  features  in  Figure  28. 

At  about  fiducial  (fid.)  no.  5760,  the  seawater  depth  was  measured  by  lead  line  to  be  12.5 
m.  The  AEM  data  indicate  about  11.7  m  of  water  under  the  ice  between  fid.  no.  5759  and 
5761 .  The  AEM  snow  and  ice  thickness  was  1 .5  m  vs  1 .62  m  of  ice  and  0.2  m  of  snow  meas¬ 
ured  by  drilling.  If  the  AEM  snow  and  ice  thickness  includes  0.1  m  of  freeboard,  then  the 
AEM  water  depth  at  fid.  no.  5760  is  about  1.4+11.7  =  13.1  m.  This  relative  depth  is  within 
10%  of  the  measured  value. 

The  airphoto-determined  topography  for  lines  3,  4  and  5  (Fig.  25)  is  shown  in  Figure  29. 

These  lines  represent  the  corridor  down  which  the  AEM  bird  was  flown.  Note  that  the  areas 
of  rough  topography  indicated  by  letters  and  numbers  in  Figure  29  agree  with  similar  peaks 
in  Figure  28.  The  airphoto  analysis  gives  an  average  snow  and  ice  elevation  for  the  three  lines 
of  5.0  m.  If  we  assume  the  floeberg  is  grounded  in  an  average  water  depth  of  1 1  m,  then  the 
ice  keel  should  also  be  this  deep.  The  average  floeberg  snow  and  ice  thickness  is  then  about 
5.0+11  =  16  m.  The  AEM  data  give  an  average  floeberg  snow  and  ice  thickness  of  12.4  m 
between  fid.  no.  5771  and  5840  and  an  average  ‘‘water  depth”  of  4.8  m.  For  reasons  previ¬ 
ously  discussed,  the  4.8  m  of  water  most  likely  represents  seawater  among  the  ice  block  keel 
structure,  which  is  in  contact  with  the  seabed.  This  is  further  indicated  by  the  very  low  con¬ 
ductivity  for  the  “seawater”  zone  under  the  floeberg  (Table  A8).  In  short,  the  ice  block  and 
seawater  conglomeration  has  affected  the  electromagnetic  response,  reducing  the  apparent 
conductivity  below  that  of  seawater  but  increasing  the  conductivity  far  above  that  of  sea  ice. 

Since  the  water  depth  under  the  floeberg  is  essentially  zero,  the  conductivity  of  the  sedi¬ 
ments  forming  the  shoal  takes  on  some  importance.  This  is  true  in  any  area  of  very  shallow 
water.  Under  this  condition  the  sediments  contribute  significantly  to  the  observed  EM  re¬ 
sponse.  If  the  water  is  deeper  than  1  skin  depth,  the  sediments  will  have  very  little  effect. 
However,  for  both  floebergs  and  first-year  pressure  ridges,  there  is  usually  some  water 
among  the  ice  blocks.  This  conductive  water  will  have  an  effect  on  the  AEM  response.  In  the 
absence  of  ground  truth  information  on  sediment  conductivity,  a  fixed  value  was  assumed, 
as  discussed  earlier. 

In  any  event,  if  the  average  4.8  m  of  “water”  under  the  floeberg  represents  a  porous  ice 
block  keel  structure,  the  average  AEM  floeberg  snow  and  ice  thickness  is  4.8  +  12.4  =  17.2 
m.  This  is  1 .2  m  thicker  than  the  relative  snow  and  ice  thickness  estimated  with  the  use  of  the 
airphoto  topographic  results  and  assumed  shoal  depth. 

Obviously  a  number  of  assumptions  were  made  to  generate  the  above  averages.  Nonethe¬ 
less,  the  exercise  is  useful  because  it  points  up  the  need  for  further  trials,  improved  data  in- 
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Figure  30.  AEM  profile  results  obtained  from  flight  6L  7.  A  portion  of  this  flight  line  is  shown  in  Figure  20. 

terpretation  techniques  and  extensive  ground  truth  measurements.  For  many  situations  the 
collection  of  ground  truth  will  not  be  a  significant  problem.  However,  for  ice  features  such 
as  the  floebergs,  detailed  ice  thickness  and  water  depth  measurements  would  be  extremely 
time  consuming  and  prohibitively  expensive.  Under  this  situation  certain  approximations, 
such  as  those  above,  will  have  to  suffice  in  assessing  the  profiling  capability  of  the  AEM  sys- 
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Flight  line  6L7  began  in  shallow  near-shore  water  and  extended  northwestward  over  the 
fast  ice  into  deeper  water.  The  location  of  a  portion  of  this  flight  track  is  shown  in  Figure  20. 
Sandpiper  and  North  Star  islands  are  man-made  gravel  exploration  islands.  Sandpiper  Island 
was  built  in  14.9  m  of  water  and  North  Star  in  13.7  m.  To  the  southeast  of  Sandpiper  Island 
the  flight  line  passed  another  man-made  structure.  Seal  Island,  built  in  1 1 .9  m  of  water.  Our 
flight  line  ran  off  the  north  side  of  the  islands,  where  the  water  was  about  O.S  m  deeper.  The 
relative  water  depths  off  Seal,  Sandpiper  and  North  Star  islands  were  about  12.S,  14  and 
15.5  m,  respectively. 

The  AEM  snow  and  ice  thickness  and  under-ice  water  depth  along  track  6L7  are  graphical¬ 
ly  presented  in  Figure  30.  The  nonpressured  snow  and  ice  thickness  is  seen  to  vary  but  aver¬ 
ages  1 .89  m  thick  between  Seal  and  Sandpiper  islands.  The  variations  in  the  snow  and  ice  re¬ 
lief  are  due  to  real  snow  and  ice  thickness  variations  as  well  as  to  system  noise  and  drift  prob¬ 
lems  and  low  laser  accuracy.  Nevertheless,  the  average  AEM  snow  and  ice  thickness  is  in  rea¬ 
sonable  agreement  with  six  drill  hole  measurements,  which  averaged  1.74  m. 

The  shoaling  of  the  seabed  at  the  location  of  the  man-made  islands  (Fig.  30)  is  apparently 
due  to  the  AEM  system  sensing  the  submerged  slope  of  the  gravel  islands.  The  average 
under-ice  water  depth  determined  by  AEM  off  Seal,  Sandpiper  and  North  Star  islands  was 
about  11.75,  13.25  and  14  m,  respectively.  Adding  the  average  AEM  snow  and  ice  thickness 
(minus  the  freeboard)  of  about  1.75  m  to  these  under-ice  water  depths  gives  relative  water 
depths  at  each  island  of  about  13.5,  15  and  15.75  m,  respectively.  Compared  to  the  relative 
values  given  above,  the  AEM  data  overestimated  the  apparent  water  depth  but  by  under 
10%.  This  will  be  further  evaluated  in  our  next  field  program. 
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CONCLUDING  REMARKS 


The  May  1985  field  program  was  undertaken  to  determine  the  feasibility  of  using  a  multi- 
frequency  airborne  electromagnetic  sounding  system  for  measuring  ice  thickness.  Funding 
and  time  constraints  necessitated  the  use  of  a  "standard”  geophysical  survey  system  and 
programs  not  fully  developed  for  data  interpretation.  Nevertheless,  the  results  proved  en¬ 
couraging  in  that  reasonable  measurements  were  made  of  both  thick  and  thin  first-year  sea 


ice  as  well  as  second-year  ice  of  varying  thicknesses.  No  other  remote  measurement  system 
has  demonstrated  the  capability  of  profiling  such  a  variety  of  ice  types  and  thicknesses.  In 
addition  to  profiling  sea  ice  thickness,  information  was  obtained,  for  the  first  time  from  a  re¬ 
mote  sensor,  on  seawater  conductivity  and  the  depth  of  water  under  a  sea  ice  cover. 

For  the  next  field  program  the  AEM  system  is  undergoing  a  major  redesign.  A  new  bird 
shell  has  been  fabricated  of  Kevlar  fabric.  This  shell  will  be  about  0.3S  m  in  diameter  and  3.5 
m  long.  The  reduced  shell  size,  along  with  new  coil  and  electronic  mounting  platforms  inside 
the  bird  and  lighter  coils,  should  reduce  the  weight  of  the  bird  to  about  half  that  of  the  bird 
used  in  1985.  The  smaller  bird  will  not  only  reduce  shipping  costs  of  the  AEM  system  to  re¬ 
mote  areas  but  more  importantly  is  a  major  first  step  to  downsize  the  bird  for  deployment 
from  a  fixed-wing  aircraft. 

The  new  bird  will  have  coils  positioned  to  better  define  ice  keels  and  their  depth.  Coil  fre¬ 
quencies  will  be  about  0.8, 4.5,  and  50  or  60  kHz.  The  0.8-kHz  data  should  allow  for  profil¬ 
ing  water  depth  to  at  least  10  m,  and  the  4.5-kHz  data  should  provide  water  conductivity  in¬ 
formation.  The  high-frequency  data  will  be  used  to  obtain  sea  ice  conductivity.  Theoretical 
modeling  indicates  that  resolution  of  sea  ice  conductivity  to  better  than  0. 1  S/m  may  be  pos¬ 
sible  when  sounding  at  50  kHz,  and  there  is  a  1-ppm  system  precision.  If  the  above  sea  ice 
conductivity  determination  proves  successful,  it  would,  in  principle,  be  possible  to  determine 
the  bulk  porosity  of  the  sea  ice  and  infer  its  strength  based  upon  the  work  of  Kovacs  et  al. 
(1987).  Reduction  of  system  noise  and  drift  are  of  paramount  importance  to  the  achievement 
of  the  desired  level  of  precision  in  sea  ice  conductivity  measurement. 

Techniques  for  reducing  system  noise  and  drift  have  been  investigated,  and  acceptable  so¬ 
lutions  will  be  implemented  in  the  new  AEM  sea  ice  thickness  measurement  system.  The  dy¬ 
namic  range  of  the  system  will  be  increased  about  32  times  over  that  in  the  1985  survey  sys¬ 
tem.  This  increase  will  be  achieved  through  the  use  of  a  16-bit  data  acquisition  system  with  a 
20- V  digitization  window.  Improvements  to  the  calibration  technique  will  be  implemented 
and  will  include  a  new  unique  built-in  absolute  phasing  and  calibration  technique.  Calibra¬ 
tion  methods  which  can  be  applied  at  the  postprocessing  phase  of  data  interpretation  have 
also  been  developed,  because  standard  airborne  electromagnetic  calibration  techniques  used 
for  mineral  exploration  proved  inadequate  for  sea  ice  thickness  measurement.  These  changes 
should  provide  the  first  truly  absolute  measurements  ever  collected  with  an  AEM  system  and 
include  the  capability  to  deal  with  negative,  missing  and  saturated  data,  as  well  as  the  ability 
to  recognize  the  best  possible  interpretation  of  the  data., 

A  Micro  VAX  II  computer  has  been  obtained  for  the  1987  field  study.  This  unit  should 
vastly  improve  the  speed  of  data  interpretation.  We  expect  a  processing  time  of  about  1  sec¬ 
ond  per  data  point.  An  array  processor  for  the  computer  is  being  considered  for  the  future, 
and  this  addition  should  increase  processing  to  at  least  20  data  points  per  second. 

The  1985  AEM  system  had  a  sampling  rate  of  10  data  points  per  second.  The  new  ice 
measurement  system  should  have  a  sampling  rate  of  20  data  points  per  second.  This  rate  pro¬ 
vides  for  an  ice  thickness  measurement  about  every  2.5  m  at  a  flight  speed  of  100  knots. 
Therefore,  with  an  array  processor  it  should  be  possible  for  the  computer  to  interpret  AEM 
sea  ice  thickness  data  in  real  time  when  the  data  are  collected  at  a  sampling  rate  of  20  per  sec¬ 
ond. 

Future  plans  include  shortening  the  AEM  bird  further  to  allow  it  to  be  deployable  by  an 
existing  fixed-wing  aircraft  launcher  system,  using  a  digital  receiver  to  replace  the  current 
analog  signal  processing  system,  developing  a  bird  that  generates  a  wideband  time-domain 
signal,  rather  than  the  discrete  continuous-wave  coils  now  used,  and  including  a  global  satel¬ 
lite  positioning  system  to  provide  highly  accurate  flight  track  location  information. 
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APPENDIX  A:  FIELD  SURVEY  RESULTS 


Table  Al.  Refrozen  lead  flight  results. 


Fiducial 

No. 

Pitch 
(Deg. ) 

AEM  Bird 
Roll 
(Deg.) 

Blev. 

(m) 

Ice 

Thiele. 

(m) 

Water 

Depth 

(n) 

Water 

Conduct 

(S/m) 

2900 

1.5 

-0.5 

34.87 

0.63 

18.48 

2.46 

2901 

1.3 

-0.6 

35.01 

0.75 

18.62 

2.47 

2902 

1.9 

-0.5 

35.29 

0.76 

18.77 

2.46 

2903 

2.1 

-0.6 

35.65 

0.71 

18.89 

2.44 

2904 

2.2 

-0.7 

35.61 

0.95 

18.95 

2.41 

2905 

2.9 

-0.1 

35.42 

1.03 

19.25 

2.41 

2906 

2.9 

-0.8 

35.40 

0.83 

19.14 

2.44 

2907 

2.6 

-0.8 

35.12 

0.80 

19.00 

2.42 

2908 

1.3 

-l.l 

34.80 

0.83 

19.18 

2.40 

2909 

1.2 

-0.4 

34.53 

1.00 

19.02 

2.41 

2910 

0.9 

-0.5 

34.75 

0.84 

19.34 

2.40 

2911 

1.2 

-0.6 

35.04 

0.71 

19.34 

2.43 

2912 

1.8 

-0.6 

34.90 

1.10 

19.09 

2.44 

2913 

1.9 

-0.9 

35.12 

1.16 

19.35 

2.43 

2914 

3.0 

-0.8 

34.97 

1.49 

19.78 

2.40 

2915 

2.4 

-0.7 

35.67 

1.13 

20.05 

2.35 

2916 

2.4 

-0.7 

36.15 

1.09 

19.72 

2.35 

2917 

2.6 

-0.9 

36.50 

1.10 

19.50 

2.37 

2918 

2.0 

-1.2 

36.78 

1.02 

18.86 

2.40 

2919 

2.3 

-0.8 

37.01 

0.77 

19.14 

2.42 

2920 

1.9 

-0.7 

36.90 

1.05 

18.46 

2.46 

2921 

2.2 

-0.7 

37.22 

0.69 

18.58 

2.46 

2922 

2.1 

-0.8 

36.95 

0.74 

18.67 

2.47 

2923 

2.0 

-0.4 

36.43 

0.67 

19.89 

2.43 

2924 

2.4 

-0.6 

35.39 

1.12 

19.04 

2.45 

2925 

2.2 

-0.6 

35.08 

0.78 

19.25 

2.43 

2926 

1.5 

-0.9 

34.58 

0.71 

19.13 

2.43 

2927 

1.1 

-0.7 

34.04 

0.63 

19.27 

2.42 

2928 

1.6 

-0.7 

33.06 

0.86 

19.06 

2.42 

2929 

1.2 

-0.7 

32.57 

0.73 

18.58 

2.45 

2930 

0.8 

-0.9 

31.93 

0.70 

18.24 

2.46 

2931 

0.8 

-0.9 

31.57 

0.65 

18.62 

2.45 

2932 

1.0 

-0.8 

31.64 

0.57 

18.90 

2.43 

2933 

1.1 

-0.5 

32.22 

0.37 

18.76 

2.42 

2934 

1.3 

-0.6 

32.76 

0.48 

18.62 

2.44 

2935 

1.0 

-0.7 

33.34 

0.65 

18.08 

2.48 

2936 

1.5 

-0.6 

33.87 

0.78 

18.63 

2.47 

2937 

1.9 

-0.4 

34.77 

0.59 

18.91 

2.44 

2938 

2.5 

-0.9 

35.39 

0.68 

19.05 

2.43 

2939 

3.3 

-0.9 

35.98 

0.66 

18.46 

2.48 

2940 

3.4 

-1.4 

36.24 

0.68 

19.09 

2.44 

2941 

3.4 

-1.0 

36.13 

0.94 

19.03 

2.44 

2942 

2.5 

-0.9 

36.19 

1. 10 

19.13 

2.45 

2943 

1.9 

-0.6 

36.75 

0.86 

19.51 

2.42 

2944 

1.5 

-0.7 

37.64 

0.44 

19.04 

2.44 

2945 

0.7 

-0.2 

37.61 

0.73 

18.73 

2.47 

2946 

-0.3 

-0.5 

37.99 

0.47 

18.55 

2.49 

2947 

0.6 

-0.5 

37.68 

0.43 

19.84 

2.43 

2948 

2.2 

-0.6 

37.13 

0.71 

19.57 

2.42 

2949 

3.0 

-0.9 

36.70 

0.99 

19.07 

2.45 

31 


m 


Table  A1  (coat’d).  Refrozen  lead  flight  results. 


rlducial 

No. 

Pitch 

(Deg.) 

AEM  Bird 
Roll 
(Deg. ) 

Elev. 

(n) 

Ice 

Thick. 

(■) 

Water 

Depth 

(■) 

Water 

Conduct 

(S/a) 

2950 

2.9 

-1.1 

36.66 

0.75 

19.19 

2.46 

2951 

3.0 

-1.0 

35.65 

1.16 

19.53 

2.45 

2952 

2.1 

-0.7 

35.39 

0.82 

19.25 

2.47 

2953 

1.4 

-0.9 

35.36 

0.60 

18.63 

2.46 

2954 

0.8 

-1.0 

35.28 

0.67 

18.29 

2.49 

2955 

0.2 

-0.4 

35.55 

0.45 

18.86 

2.44 

2956 

0.0 

-0.6 

35.92 

0.35 

18.63 

2.47 

2957 

0.4 

-0.7 

35.84 

0.62 

18.96 

2.45 

2958 

1.9 

-0.5 

36.06 

0.57 

19.18 

2.45 

2959 

2.2 

-0.9 

35.86 

1.11 

18.68 

2.47 

2960 

2.8 

-0.7 

36.31 

0.86 

19.09 

2.46 

2961 

2.4 

-0.9 

36.23 

1.04 

19.28 

2.44 

2962 

2.1 

-0.7 

36.28 

0.91 

19.14 

2.45 

2963 

2.1 

-0.9 

35.95 

1.13 

18.59 

2.47 

2964 

2.0 

-0.9 

36.18 

0.82 

18.41 

2.49 

2965 

1.3 

-0.7 

36.50 

0.65 

18.83 

2.45 

2966 

0.1 

-0.7 

36.76 

0.86 

18.83 

2.46 

2967 

0.3 

-0.2 

37.52 

0.66 

19.03 

2.45 

2968 

1.3 

-0.6 

38.37 

0.58 

18.64 

2.48 

2969 

1.4 

-0.7 

38.60 

0.99 

19.02 

2.45 

2970 

1.9 

-0.5 

39.18 

0.75 

19.25 

2.47 

2971 

2.6 

-0.7 

38.88 

0.88 

19.77 

2.46 

2972 

3.1 

-0.8 

38.58 

1.02 

19.21 

2.48 

2973 

3.0 

-1.1 

38.49 

0.81 

19.35 

2.46 

2974 

2.3 

-0.9 

37.73 

0.90 

19.53 

2.44 

2975 

1.5 

-0.7 

37.46 

0.62 

19.10 

2.45 

2976 

1.3 

-0.6 

36.81 

0.78 

18.56 

2.48 

2977 

-0.3 

-0.8 

36.50 

0.65 

19.04 

2.47 

2978 

0.6 

-0.4 

36.25 

0.65 

19.04 

2.47 

2979 

1.2 

-0.6 

35.83 

0.95 

18.74 

2.47 

2980 

1.1 

-0.9 

35.87 

0.80 

18.92 

2.47 

2981 

1.9 

-0.5 

35.76 

0.80 

18.92 

2.47 

2982 

1.9 

-0.7 

35.68 

0.75 

19.38 

2.45 

2983 

1.9 

-1.0 

35.35 

0.83 

19.01 

2.46 

2984 

1.9 

-1.0 

34.90 

0.84 

19.10 

2.46 

2985 

1.5 

-0.7 

34.41 

0.84 

19.15 

2.45 

2986 

0.7 

-0.5 

34.12 

0.78 

18.90 

2.46 

2987 

0.5 

-0.7 

34.25 

0.47 

18.77 

2.46 

2988 

0.4 

-0.7 

33.76 

0.81 

18.74 

2.45 

2989 

1.0 

-0.7 

33.85 

0.80 

18.65 

2.46 

2990 

1.2 

-0.5 

34.06 

0.68 

18.90 

2.46 

2991 

1.5 

-0.6 

34.04 

0.81 

18.74 

2.47 

2992 

2.1 

-0.7 

34.04 

0.86 

19.04 

2.46 

2993 

2.6 

-0.7 

34.10 

0.85 

19.32 

2.44 

2994 

2.5 

-0.7 

34.18 

0.90 

19.15 

2.45 

2995 

2.3 

-0.8 

34.44 

0.80 

16.92 

2.47 

2996 

2.0 

-0.6 

34.63 

0.68 

19.25 

2.44 

2997 

1.7 

-0.6 

34.78 

0.72 

18.99 

2.46 

2998 

1.9 

-0.9 

34.91 

0.90 

18.66 

2.47 

2999 

1.9 

-0.7 

35.29 

0.76 

18.83 

2.46 

3000 

2.5 

-0.8 

35.46 

0.84 

18.82 

2.47 

>v 

£ 

^  m 
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Table  A2.  Survey  line  station  driU-hoie-measored  snow  and  ice  thicknesses. 


Total 

Snow  Ice  Snow  and  Ice 

Depth  Thickness  Thickness 
(cm) _ (■) _ (n) 

a.  Survey  line  A  station. 


0 

7 

1.78 

1 

20 

1.80 

2 

42 

1.28 

3 

27 

1.58 

4 

20 

2.25 

5 

8 

2.97 

6 

36 

2.12 

7 

40 

2.50 

8 

17 

2.31 

9 

25 

2.75 

10 

10 

3.50 

11 

38 

3.62 

12 

8 

4.37 

13 

7 

3.93 

14 

28 

4.32 

14.5 

40 

4.80 

15 

50 

4.20 

15.5 

46 

4.59 

16 

10 

5.95 

16.5 

5 

7.65 

17 

29 

6.61 

17.5 

31 

6.29 

18 

21 

5.59 

19 

9 

4.51 

20 

28 

3.07 

21 

9 

3.51 

22 

8 

3.67 

23 

47 

3.13 

24 

41 

2.99 

25 

41 

2.99 

26 

5 

3.25 

27 

23 

3.07 

28 

8 

3.82 

29 

22 

2.43 

30 

5 

1.95 

31 

24 

2.06 

32 

9 

2.81 

b.  Survey  line  B  station. 


Total 

Snow  Ice  Snow  and  Ice 

Depth  Thickness  Thickness 

(cm) _ (■) _ (n) _ 


c.  Survey  line  C  station. 

20  3.40 

8  3.19 

10  1.60 


3 

o 

9 

1.91 

2.00 

17.5 

15 

6.05 

i 

10 

2.50 

2.60 

18 

21 

4.84 

V 

2 

12 

4.43 

4.55 

19 

24 

3.76 

t# 

3 

48 

3.02 

3.50 

20 

55 

3.30 

4 

18 

2.27 

2.45 

21 

6 

3.79 

5 

4 

2.06 

2.10 

22 

34 

2.76 

6 

39 

2.71 

3.10 

23 

11 

3.19 

7 

18 

3.52 

3.70 

24 

20 

2.55 

!  ■ 

8 

7 

3.63 

3.70 

25 

39 

2.31 

«j 

9 

7 

3.58 

3.65 

26 

8 

3.42 

*( 

10 

22 

2.98 

3.20 

27 

43 

2.77 

r 

11 

20 

3.40 

3.60 

28 

9 

3.71 

* 

12 

9 

4.11 

4.20 

29 

9 

4.16 

13 

10 

4.10 

4.20 

30 

9 

3.41 

1 

14 

6 

4.64 

4.70 

31 

15 

2.85 

i 

32 

38 

2.72 

Table  A4.  Line  6L4  flight  results 


Tabic  A5.  Line  8L2  flight  results, 


AEM  Bird 

Snow-Ice 

Hater 

Fiducial 

Pitch 

Roll 

Elev. 

Thickness 

Conduct 

No. 

(Deg.) 

(Deg.) 

.  (■) 

(■) 

(S/e) 

2497 

5.7 

-0.3 

45.76 

2.93 

2.52 

2498 

5.9 

-0.7 

45.99 

2.58 

2.39 

2499 

3.8 

-1.0 

44.13 

2.90 

2.27 

2500 

3.1 

-1.5 

40.87 

4.67 

2.29 

2501 

3.4 

-1.5 

39.18 

4.13 

2.44 

2502 

4.0 

-1.0 

40.52 

3.04 

2.45 

2503 

5.0 

-0.3 

42.12 

2.47 

2.45 

2504 

5.3 

0.2 

42.29 

2.71 

2.44 

2505 

4.9 

0.3 

40.56 

3.88 

2.46 

2506 

4.2 

0.0 

38.96 

4.55 

2.54 

2507 

4.9 

0.2 

38.48 

4.77 

2.63 

2508 

6.0 

0.1 

39.27 

4.20 

2.58 

2509 

7.2 

-0.5 

40.01 

4.09 

2.54 

2510 

7.0 

-0.5 

40.08 

4.18 

2.52 

2511 

5.9 

-0.7 

39.23 

4.74 

2.57 

2512 

5.4 

-1.2 

38.48 

5.15 

2.60 

2513 

5.8 

-1.1 

39.45 

4.81 

2.55 

2514 

6.3 

-0.9 

41.89 

4.01 

2.45 

2515 

6.1 

-0.2 

43.98 

3.24 

2.41 

2516 

4.8 

0.0 

43.67 

3.83 

2.51 

2517 

4.5 

0.1 

42.20 

4.90 

2.66 

2518 

5.0 

0.2 

41.76 

5.55 

2.87 

2519 

6.1 

-0.1 

43.54 

5.06 

2.98 

2520 

7.7 

-0.5 

45.99 

3.44 

2.76 

2521 

7.9 

-0.8 

45.79 

3.77 

2.72 

2522 

5.6 

-0.9 

44.00 

4.62 

2.63 

2523 

4.7 

-1.3 

41.95 

5.15 

2.63 

2524 

4.7 

-1.3 

41.62 

4.94 

2.66 

2525 

5.1 

-0.9 

43.96 

3.74 

2.64 

2526 

4.5 

0.0 

46.56 

2.79 

2.67 

2527 

4.0 

0.2 

45.96 

3.73 

2.71 

2528 

4.0 

0.2 

43.87 

4.73 

2.72 

2529 

4.6 

-0.1 

43.48 

4.97 

2.87 

2530 

5.5 

0.0 

44.19 

5.13 

2.87 
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Table  A6.  Line  8L3  flight  results, 


Fiducial 

No. 

Pitch 

(Deg.) 

AEM  Bird 
Roll 
(Deg.) 

Elev. 

(m) 

Snow-ice 

Thickness 

(a) 

Water 

Conduct 

(S/m) 

2639 

5.2 

-0.7 

37.55 

2.78 

2.46 

2640 

5.7 

-0.8 

37.93 

2.77 

2.45 

2641 

5.7 

-0.9 

38.68 

2.65 

2.43 

2642 

5.7 

-0.4 

39.04 

2.60 

2.38 

2643 

5.7 

-0.3 

38.69 

2.70 

2.30 

2644 

5.6 

0.1 

37.67 

3.00 

2.23 

2645 

5.7 

-0.1 

37.00 

3.13 

2.26 

2646 

5.3 

-0.1 

37.34 

2.78 

2.24 

2647 

5.1 

-0.4 

37.31 

3.23 

2.39 

2648 

4.4 

-0.7 

37.12 

3.52 

2.44 

2649 

3.4 

-0.8 

36.83 

3.34 

2.33 

2650 

2.8 

-0.6 

36.05 

3.57 

2.23 

2651 

3.8 

-0.9 

35.52 

3.95 

2.14 

2652 

5.3 

-0.8 

35.69 

4.08 

2.10 

2653 

6.5 

-0.7 

37.24 

2.96 

2.03 

2654 

5.9 

-0.3 

37.51 

3.28 

2.05 

2655 

5.4 

0.2 

37.66 

2.96 

2.10 

2656 

4.8 

0.2 

36.26 

4.02 

2.11 

2657 

4.4 

0.0 

36.85 

3.26 

2.16 

2658 

4.0 

-0.6 

37.20 

2.83 

2.22 

2659 

3.8 

-0.8 

36.91 

2.80 

2.18 

2660 

3.4 

-l.l 

35.93 

3.12 

2.21 

2661 

2.9 

-0.9 

35.09 

3.26 

2.23 

2662 

3.8 

-0.7 

35.44 

2.65 

2.19 

2663 

5.1 

-0.7 

36.28 

2.26 

2.13 

2664 

5.5 

-0.4 

36.74 

2.23 

2.11 
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Table  A7.  Line  8L4  flight  results. 


AEM  Bird 


Snow-ice  Water 


*v: 

K'f.- 


V 


v\»  v* 


Table  A9.  Floeberg  C  flight  results 


Table  A9  (coat’d).  Floeberg  C  flight  malts, 


AEM  Bird 

Ice 

Water 

Water 

Fiducial 

Pitch 

Roll 

Elev. 

Thick. 

Depth 

Conduct 

No. 

(Deg.) 

(Deg.) 

(a) 

(a) 

(a) 

(S/a) 

5796 

4.4 

-1.4 

31.22 

16.10 

3.42 

1.27 

5797 

5.0 

-1.7 

34.64 

14.35 

4.02 

1.20 

5798 

5.7 

-1.1 

35.54 

14.35 

4.02 

1.20 

5799 

5.8 

-0.7 

36.90 

12.81 

3.91 

1.20 

5800 

4.9 

-0.4 

37.47 

12.81 

3.91 

1.20 

5801 

4.2 

-0.4 

34.30 

13.64 

3.73 

1.20 

5802 

4.3 

-0.4 

33.37 

13.80 

3.86 

1.20 

5803 

5.3 

-0.9 

36.02 

11.67 

3.64 

1.20 

5804 

5.8 

-0.9 

39.56 

9.43 

3.10 

1.20 

5805 

5.9 

-1.1 

39.68 

9.43 

3.10 

1.20 

5806 

5.5 

-1.3 

36.82 

12.10 

3.87 

1.00 

5807 

4.9 

-1.6 

33.21 

15.18 

3.93 

1.37 

5808 

5.7 

-1.8 

34.89 

13.92 

3.67 

1.45 

5809 

6.0 

-1.4 

37.15 

13.92 

3.67 

1.30 

5810 

5.8 

-1.0 

39.35 

12.49 

3.61 

1.32 

5811 

4.8 

-0.4 

38.67 

13.28 

3.73 

1.38 

5812 

3.9 

-0.3 

38.96 

13.28 

3.73 

1.38 

5813 

4.4 

-0.5 

38.14 

12.42 

3.87 

1.38 

5814 

5.2 

-0.6 

40.23 

12.42 

3.87 

1.38 

5815 

6.6 

-1.2 

46.46 

6.87 

4.24 

1.38 

5816 

6.9 

-1.6 

48.19 

5.55 

4.23 

1.38 

5817 

5.8 

-1.5 

45.12 

8.88 

4.22 

1.38 

5818 

5.0 

-1.6 

43.04 

10.51 

4.01 

1.38 

5819 

5.2 

-1.7 

41.48 

12.44 

3.60 

1.38 

5820 

5.6 

-1.5 

43.59 

12.44 

3.60 

1.38 

5821 

4.9 

-0.8 

45.99 

10.49 

3.66 

1.38 

5822 

4.1 

-0.1 

44.23 

10.49 

3.66 

1.38 

5823 

3.2 

-0.1 

42.11 

11.22 

3.83 

1.38 

5824 

3.7 

-0.3 

38.02 

12.77 

3.74 

1.38 

5825 

4.3 

-0.8 

36.18 

13.25 

5.25 

1.38 

5826 

5.6 

-1.0 

38.26 

11.03 

4.79 

1.38 

5827 

6.5 

-1.5 

40.36 

8.43 

4.69 

1.38 

5828 

6.0 

-1.4 

40.28 

7.83 

4.96 

1.38 

5829 

5.3 

-1.5 

37.04 

9.97 

5.29 

1.38 

5830 

4.8 

-1.8 

37.94 

7.86 

6.17 

1.38 

5831 

4.7 

-1.6 

38.63 

6.44 

6.30 

1.38 

5832 

4.7 

-1.0 

40.40 

4.74 

6.34 

1.38 

5833 

4.6 

-0.7 

41.07 

3.64 

6.89 

1.38 

5834 

4.1 

-0.1 

39.57 

3.64 

6.89 

1.14 

5835 

3.8 

-0.3 

37.54 

3.70 

7.25 

1.28 

5836 

4.3 

-0.4 

35.49 

3.73 

7.46 

1.40 

5837 

5.7 

-0.8 

35.52 

2.95 

7.77 

1.45 

5838 

6.4 

-1.2 

35.37 

3.10 

8.51 

1.55 

5839 

6.4 

-1.3 

35.11 

2.90 

8.85 

1.60 

5840 

5.4 

-1.3 

33.71 

3.04 

9.09 

1.67 
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A  facsimile  catalog  card  in  Library  of  Congress  MARC  format  is  repro¬ 
duced  below. 


Kovacs,  Austin 

Airborne  electromagnetic  sounding  of  sea  ice  thickness  and  sub-ice  bathy¬ 
metry  /  by  Austin  Kovacs,  Nicholas  Valleau  and  J.  Scott  Holladay.  Han¬ 
over,  N.H.:  U.S.  Army  Cold  Regions  Research  and  Engineering  Labora¬ 
tory;  Springfield,  Va.:  available  from  National  Technical  Information  Ser¬ 
vice,  1987. 

iv,  48  p.,  illus. ;  28  cm.  (CRREL  Report  87-23.) 

Bibliography:  p.  29. 

1.  Airborne  electromagnetic  sounding.  2.  Ice  thickness.  3.  Remote  sens¬ 
ing.  4.  Sea  ice.  I.  I.  Valleau,  Nicholas.  11.  Holladay,  J.  Scott.  III.  United 
States  Army.  Corps  of  Engineers.  IV.  Cold  Regions  Research  and  Engineer¬ 
ing  Laboratory.  V.  Series:  CRREL  Report  87-23. 
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